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A new type of hidden order in many body systems is explored. This order appears in states which 
are analogues to charge density waves, or spin density waves, but involve anomalous particle-hole 
correlations that are odd in relative time and frequency. These states are shown to be inherently 
different from the usual states of density waves. We discuss two methods to experimentally observe 
the new type of pairing where a clear distinction between odd and even correlations can be detected: 

(i) by measuring the density-density correlation, both in time and space and (ii) via the conduc¬ 
tivity which, according to the Kubo formula, is given by the current-current correlation. An order 
parameter for these states is defined and calculated for a simple model, illuminating the physical 
nature of this order. 


I. INTRODUCTION 

The conventional way to consider ordered many body 
systems is the Landau-Ginzburg formalism of phase tran¬ 
sition that is centered on order parameter as the attribute 
of order. That appearing order is classified as the anoma¬ 
lous correlations that describe the emergent state. For 
example, in charge density wave the order parameter is 


wave correlations. This odd time density waves would 
describe dynamic order that does not manifest itself in 
any static density wave. The situation can be viewed 
as a dynamic order that is hidden from the conventional 
spectroscopies of charge and spin density waves. 


II. INTRODUCING ODD TIME CDW 


Aq(k) — (fl^q^Ok.a)' (1) 

where aj^ is creating (annihilating) an electron with mo¬ 
mentum k and spin a. Similarly, the order parameter for 
spin density waves is given by 

A q (k) = (oj c _|_ qa Cl] i) ^)iT Q , | g. (2) 

where a is a vector of Pauli matrices. Summation over 
spin indices is implied, both in 0 and 0 and these 
anomalous correlations are taken at equal time for the 
involved operators. Static values as written are often 
taken as the order parameters of the charge (spin) density 
waves 

As our view of correlations evolves, it becomes ap¬ 
parent that materials can exhibit composite meta or¬ 
der that significantly expands the old Ginzburg-Landau 
paradigm. Examples of new orders include topological 
order with no local order parameter®®. Another ex¬ 
tension of the concept of order is the notion of odd 
frequency or odd-time correlations, first proposed by 
Berezinskii.' 8 Odd-frequency order in superconductors 
and superfluids has been proposed for numerous real¬ 
izations and most likely to occur i n su perconducting 
heterostructures, in nanoscale device^® and in mulit- 
band superconductortPd. In addition to superconductiv¬ 
ity odd time orders were expanded to spin nematiciP 2 and 
BECP. Here we wish to further advance the notion of 
odd time orders by considering the odd in time density 


The odd time, or odd frequency, orders are character¬ 
ized by correlations that vanish when (a^a) expectaion is 
taken at equal time. We now focus on charge and spin 
density waves (CDW,SDW). Since the discussion, in our 
context, for these two types of states, is highly similar, 
for the sake of simplicity, we will first treat only CDW, 
described by eq. ([l]). Since spin indices are summed over, 
they are suppressed in all expressions. In order to treat 
SDW, described by eq. 0 one simply has to reintro¬ 
duce the spin indices and multiply by a throughout the 
derivation. 

We can generalize the correlation |l]) to include oper¬ 
ators ak(7") which act in Matsubara time r. The result 
is the anomalous green function 


A q( k U) 



( 3 ) 


where T t is the time ordering operator and 9(x) is the 
Heaviside step function. Since the system is assumed to 
be invariant to time translation, A q (k, r) is a function of 
the relative time only. 

There can be a few possibilities for the time depen¬ 
dency of A q (k, t) at short times, r —> 0. If there is a 
zero time order, i.e. a time independent part, then we 
recover the usual form of CDW with non vanishing order 
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parameter given by |T]). Then the time dependence of 
the correlator at nonzero times would reflect the retar¬ 
dation effects that occur at a later time. If, on the other 
hand, the equal time part vanishes, but the first order 
A q (k, r) ~ t remains finite, the situation is different. In 
such a case, all the well-known results regarding CDW 
might not be valid, but nonetheless there will be non¬ 
trivial density-density correlations. We are thus facing 
the possibility to have a nontrivial correlations that have 
no equal time signature, i.e. a hidden order. 

To capture the dynamic nature of the correlation we 
consider a scenario is when A q (k, r) is odd in time: 

A q (k, t) = —A q (k, — t). (4) 

In this case it is clear that A q (k) = A q (k, 0) = 0. This 
implies that if the system is described by quantity obey¬ 
ing a condition of the form Q, it is in a state that pos¬ 
sess a dynamic kind of order. This state might share 
many of the characteristics, of the usual CDW, by it is 
fundamentally different in that there is no static density 
correlations. 


III. DETECTION OF ODD TIME CDW 

Here we outline experimental observation that would 
test the odd frequency density order. The usual type 
of CDW, which is described by (|T|) is manifested with a 
modulation of the density 

n(r) = ('I' t (r)'k(r)) = n 0 + n q cos(qr + </>), (5) 

where T(r) = J dke !kr ak, no is the average density, n q 
is the amplitude of the modulation with wave number q 
and (f> is an arbitrary phase. When writing Eq. © we 
ignore any density modulations which are not directly 
related to CDW, so we explicitly assume a uniform den¬ 
sity for n q = 0. Indeed, the uniform density n(r) = no 
would be the case when the equal time order parameter 
vanish, A q (k) = 0, as would occur, for example, in the 
odd frequency case. So one has to study higher order cor¬ 
relations to reveal this kind of order. This can be done 
by looking at the density correlation, both in space and 
time 

(T T h(r, r)h(r', 0)) = n 2 (6) 

+ (T T ^(r, (r', 0 )><T t * V, 0)tf (r, r)), 

where we used Wick’s theorem and assumed that n = 
(n(r, t)) is spatially and temporally independent. This 
assumption involves two important aspects. First the 
spatial homogeneity means we do not have the usual 
CDW state, which is manifested by the modulations 
([5]). Second the time r in (n(r,r)) is the absolute 
time, or “center of mass” time, as opposed to r in 
(T T n(r, r)n(r , ,0)} which represent a relative time, or 
some correlation time. Dependency on the relative time 
can comes from microscopical dynamics. Here we note 



Figure 1. (Color online) An illustration of the typical sit¬ 
uation of CDW. Correlations between electrons on opposite 
sides of the Fermi surface appear and a gap is opened, around 
the Fermi energy, proportional to the correlation magnitude. 
Since the involved electrons should be in the vicinity of the 
Fermi wave vector ±fc.F, the momentum difference, which also 
determine the wave vector of the density modulations, have 
to be jq| = 2k,F■ The correlation operator can be written 
either as annihilation of electron at fop and creation at — Icf 
so A q (k) = {a^_ kF ak F ), which make q = — 2fcp, or vice versa 
Aq(k) = (a^a-fcp), which makes q = 2ftp. Note that for 
odd time CDW A q (k) = A q (k, r = 0) = 0 and a gap might 
not open at fcp, but it is still reasonable to assume that cor¬ 
relations with the same momentum properties appear. This 
assumption is represented in eq (|8|. 


that there are discussions on states that depend on the 
overall time. Such a scenario can imply the system is 
not in a st eady state, or even the formation of a time 
crystafSE^. We do not address those issues in this work. 

The quantity of interest for us is the second term on the 
right hand side of ([g]) , x(r, r',r) = (T r ?f(r, r)n(r', 0)) — 
n 2 . By going over to momentum space we can write it 
using (|3]) : 

X(r,r , ,r)= J A q (k, r)A q ,(k', -r), 

(?) 

where r = r — i - ', k = k — k' and f dIC = dkdk'dqdq' 
means integration on all momentum variables. Let us 
now consider a case similar to the typical scenario for 
CDW, shown in Fig [I] where the abnormal correlations 
occur on the two opposite sides of a Fermi surface. This 
implies q = —2k = ±2k/, where k/ is the Fermi wave 
vector. So the order parameter is given by 

A q (k, r) = A+(r)5(k - k/)c5(q + 2k/) 

+ A_(T)5(k + k / )*(q-2k / ). (8) 

where A + (r) and A_ (r) are in principle independent 
quantities. Here, for simplicity we choose d-function as 
the form of A q (k, t). A more realistic assumption would 
be that it is sharply peaked function at those values, so 
that A(t) = J dkdqA q (k, r) with integration inter¬ 
vals are around the points k = (—)k/, —q = (—)2k/. 
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Inserting ^ in Q we have, 

X(r,r',r) = A +(r) [A + (-r)e s4Rk/ + A_(-r)] 

+ A_(r) [A_(—r)e _i4Rk/ + A+(-r)] . (9) 

The fact that x(r, r', r) depends on the “center of mass” 
coordinate R = (r + r')/2 is a clear manifestation of 
the broken translation symmetry associated with CDW. 
Since we are considering a situation where the average 
density is uniform, the broken symmetry that can be ob¬ 
served in the density-density correlation is a strong indi¬ 
cation to a hidden order residing in the system. 

The first two terms in ([TJ) are oscillating with effec¬ 
tive wave vector of 4 kf. Such terms would appear for 
the usual case of CDW, if one looks at density-density 
correlations, in addition to terms with a wave vector of 
2 kf, which is also how the density is modulated. The 
more rapid oscillations of might be difficult to observe in 
certain systems. Neglecting the oscillating term we are 
left with a constant, i.e. r independent term. In the sim¬ 
ple case where A+(t) = A_(r) = A(r), the sign of this 
term is given by the time parity of A(t), regardless of the 
sign of A(r) itself. Thus, if A(t) = —A(—r), there will 
be a negative contribution from x(r, r 7 , t) to the density- 
density correlation. We can consider a scenario, where 
this fact can be the basis for experimental evidence for 
odd time correlation, when the density-density correla¬ 
tion is a measured function of some experimental param¬ 
eter, like temperature. If the correlation is constant in 
some region and then starts changing, the direction of 
the change can tell us whether A q (k, r) is odd or even in 
time, i.e. a decrease in correlation would suggest an odd 
time CDW. 

Another quantity that can be of interest is the current- 
current correlations, which according to the Kubo 
formula 16 can yield the conductivity. Assuming the av¬ 
erage current vanishes (j) = 0, the current-current cor¬ 
relations are given by (see appendix | for details): 


(T T ji (r, t )jj (r', 0)) = J d/Ce _ ^ fk+rq+r q ^ x 

(k + k' + q)j(k + k' + q')j- 


Am 2 


A q (k, r)A q /(k', —r). 

( 10 ) 


Following the same procedure we used above for 
X(r, r',r), we consider the typical case, shown in F ig |T| 
where A q (k,r) is given by (|8|. Inserting ^ into (101, 
we get 


O'MMM r',0)) = 

(Mi (Mj 


[A+(t)A_(-t) + A_(t)A + (- 7 


( 11 ) 


It is interesting that the rapidly oscillating terms, which 
appear in eq 0 . vanish here because they correspond to 
zero momentum. The disappearance of a conductivity 


P T 

A+ A 

+ 
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1 
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Table I. The values of the density-density and the current- 
current correlation for di ffere nt parity behaviors of A q (k, t) 
according te eq. jf)]) and (111. A 6 denotes Aq,(/3t), P and T 
are the reversal operator for time and momentum respectively. 
The last line refer to a situation where P and T do not have 
a defined value each, but are interdependent, which is noted 
as Entang(led). From the observation of these values can one 
can learn the type, and especially the parity, of electron-hole 
correlations present in a system. 


term, accompanied by the appearance of oscillating den¬ 
sity terms seems similar, and might be closely related, to 
the usual scenario of CDW where the conductivity vanish 
together with the appearance of density modulation. 

In the same way we consider A q (k, r) to be an odd 
or even function of time, we can also consider whether 
it is odd or even in momentum. For the case de¬ 
scribed by (18]), even (odd) parity in momentum implies 
A + (r) = (A)A_(r). These properties of A q (k,r) can 
be formulated by considering A q (k, r) to be a eigenfunc¬ 
tion of two operators: the time reversal T, defined as 
TA q (k, t) = A q (k, —r) and the momentum reversal, de¬ 
fined as as PA q (k,r) = A_ q (—k, r). For each oper¬ 
ator, an eigenvalues of -1 or 1 implies the function is 
odd or even respectively. The values of x( r > r \ T ) and 
(j(r , r)j(r', 0)), for the different parity options are shown 
in Table |U 

Apart from the four parity options for the operators P 
and T, there can be another kind of behavior. For ex¬ 
ample, if A + (r) is even in time and A_(r) is odd. Such 
a situation, where the parity in one variable depends on 
the sign of the other variable suggest that the time and 
momentum parities are entangled (not necessarily in the 
quantum mechanics sense). Thus the parity of each indi¬ 
vidual operator is not defined and we refer to this state 
as entangled. There could be several possibilities for such 
a state and the values of x(r,r',T) and (j(r,r)j(r',0)), 
for one of those possibilities are shown in the last line of 
Table U 

The typical case of CDW require that A+(t) = A_(r), 
i.e. P even, and independent of the sign of the mo¬ 
mentum. One can allow to have explicit non-trivial de¬ 
pendence on the sign of the momentum i.e. A + (r) = 
—A_(t), which means we are considering p-wave CDW. 
Hence, we have both possibilities and Table [T] represent 
a general symmetry classification. 
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IV. DEFINITION OF AN ORDER PARAMETER 
FOR ODD TIME CDW 


electrons, as given by ([5]), and write an expression for n q 
in term of averages on creation operators 


If indeed this state represent a new phase, which pos¬ 
sess some order, it is described by an order parameter. 
Since the order parameter, as it is most commonly de¬ 
fined, vanish in the case of odd time correlations, we 
define a new order parameter. This can be done, for 
example, by considering the time derivative of the corre¬ 
lation at equal time, i.e. r = 0. Let us define our order 
parameter aJl 7 ^ 

D q (k) = ^A q (k,r)| T=0 . (12) 

While D q (k) does not depend on time it capture some of 
the dynamics of the state. For any model Hamiltonian 
H, we can use the Heisenberg equation ■£-0(t ) = [H, 0\ 
to write this order parameter as 


n a = e 


_ pi<t> 


dk ( a k+q fl k) = / dkA q (k). (16) 


The presence of modulations n q 7 ^ 0, implies a non 
vanishing order parameter regardless of the microscopic 
model but from ( p~5] ) it is clear that an ( 6 q ) ~ A q (k). 
So the traditional formalism can be regarded as a pos¬ 
sible recipe to construct a microscopic model which can 
explain our phenomenological quantities. 


D q (k). We have 


We can use (15) for getting an explicit expression for 


-Dq(k) — /^^5 q ',k+q (&l q ' + kq') °k+q+q' °k \ • (1") 


£> q (k) = ([ff,4 +q ]a k ). (13) 

In case the part that is most significant is the correlations 
at close to equal time r <C 1 , it is plausible to approxi¬ 
mate this correlations to first order as 


A q (k, r) ~ D q (k)r 


(14) 


and then we can learn from the order parameter (13) 


relevant correlations that are responsible for odd time 
order (|3|. 


It can be very instructive to calculate (13) for a spe¬ 


cific model and a good candidate would be the textbook 
model for CDW. However, our discussion of CDW, which 
involved a phenomenological approach and was focused 
on electron correlation, is somewhat different from the 
traditional treatment of CDW, which typically includes 
a microscopic model and is focused on average of phonon 
operators. Thus we briefly review the traditional formal¬ 
ism and show how it connects to the one we used. 

The traditional derivation of CDW^ 

involves the 

electron-phonon interaction 


Hi 


qF [bl 


k,q 


+ 6q) 


a k+q« k 


(15) 


where 6 q * is creating (annihilating) a boson with momen¬ 
tum q and g q k is a coupling function. Using mean field 
methods, a self consistent equation can be obtained for 
( 6 q ) with a solution having ( 6 q ) 0 , which represent a 

periodic distortion in the underlying lattice. The average 
(bq) is related to the gap that is opened in the electron 
spectrum and it is interpreted as an order parameter. 
From a phenomenological point of view, if one is inter¬ 
ested in the electronic structure, it is not necessary to 
consider the interaction (15), or even to involve phonons, 


in order to derive an order parameter. One can simply 
consider periodic spatial modulation n q in the density of 


In the typical case, which we discussed above, q = q' = 
2 kf. Thus the order parameter is given by electron-hole 
correlation with a momentum difference of 4 kf. This ob¬ 
servation is consistent with the behavior of y(r, r', r) and 
(j(r, r)j( r ', 0)) calculated above. Indeed this suggests 
that the hidden order we are interested in is manifested 
in higher order correlation between electron and lattice 
degrees of freedom. 


V. CONCLUSION 


In conclusion, we present the odd time (frequency) 
CDW and SDW as an extension to the notion of den¬ 
sity waves considered to date. Main feature of these odd 
time density waves is that there are no equal time modu¬ 
lation in the expectation values of spin and charge densi¬ 
ties. Yet as is clear from the overall structure, there are 
nontrivial dynamic correlations. These nontrivial corre¬ 
lations, Eq. |9]) and (11) could serve as a test for the 
odd time density wave. As such this work represents the 
extension of odd-frequency state classification to density 
waves, beyond superconductors and BEC. 


It remains to be seen what are the realistic fermion 
interactions that would enable the odd time density cor¬ 
relations. We also point out that these dynamic correla¬ 
tions with no equal time expectation value do correspond 
in a way to a “hidden order” that cannot be observed as 
the traditional charge/spin modulation. As such it can be 
a candidate for materials with hidden orders where one 
observes thermodynamic features of a transition while 
no order parameter consistent with the transition can be 
inferred from the current measurements. 
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Appendix: Current-current correlation 


The current operator is given by 

ji( r ) = (5' t (r)V J ;5'(r) - ^'(r)V i ^ t (r)) 

= J dkdk'e^'-Vaiav^^. 

The current-current correlation is given by 
(T T ji (r, t )jj (r 1 2 3 4 5 6 7 8 9 10 ,0)} = 

[ ^ c ^ r l k l~ k 2) c » r '( k 3~k4) (^1 

J 2m 2m 

(Tral ( r)a k2 (r)aj C3 ( 0 )a kj ( 0 )). 


Using Wicks theorem we can write the average on four 
operators as the product of two averages on two opera¬ 
tors. Since we are not considering the situation of Cooper 
pairs, averages containing two annihilation operators or 
two creation operators will vanish. 

Thus we are left with two ways of pairing 
the operators: (: T T a\. (V)ak 2 (r))(T r aj c (0)ak 4 (0)) and 
(A.l) (^Va ki (r)ak 4 (0))(T T aj C 3 (0)a k2 (r)). . 

The first way of pairing will simply give us 
(j( r, r))(j(r', 0 )), and we can assume for simplicity that it 
vanish. The second way of pairing is related to A q (k, r). 
We can make this relation concrete with the mapping 



k = k 4 , 

X 

►Q 

II 

(A.2) 

k' = k 2 , 
q' = k 3 - k 2 


where f dlC = dkidk 2 dk 3 dk 4 means integration on all 
momentum variables. 


Rewriting the current-current correlation above using 
these variables yields eq ( 10 ) in the main text. 
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